arXiv:hep-ph/9808203vl 2 Aug 1998 


IFT-P.049/98 

hep-ph/9808203 


Update on 


Atmospheric Neutrino^* 


M. C. Gonzalez-Garcia 


Instituto de Fisica Corpuscular - C.S.I.C/Universitat de Valencia. Spain 

and 

Instituto de Fisica Teorica, Universidade Estadual Paulista 
Rua Pamplona 145,01405-900 Sdo Paulo, Brazil 
E-mail: concha.gonzalez@uv.es 


ABSTRACT 

I summarize here the results of a global fit to the full data set corresponding 
to 535 days of data of the Super-Kamiokande experiment as well as to all other 
experiments in order to compare the two most likely solutions to the atmospheric 
neutrino anomaly in terms of oscillations in the Vt and —> Vg channels. 

1. Introduction 

Atmospheric showers are initiated when primary cosmic rays hit the Earth’s atmo¬ 
sphere. Secondary mesons produced in this collision, mostly pions and kaons, decay 
and give rise to electron and muon neutrino and anti-neutrinos fluxes H. There has 
been a long-standing anomaly between the predicted and observed /z/g ratio of 
the atmospheric neutrino fluxes i. Although the absolute individual or z/g fluxes 
are only known to within 30% accuracy, different authors agree that the z/^ /z/g ratio 
is accurate up to a 5% precision. In this resides our confidence on the atmospheric 
neutrino anomaly (ANA), now strengthened by the high statistics sample collected at 
the Super-Kamiokande experiment □. The most likely solution of the ANA involves 
neutrino oscillations. In principle we can invoke various neutrino oscillation channels, 
involving the conversion of z/^ into either z/g or (active-active transitions) or the 
oscillation of z/^ into a sterile neutrino Vg (active-sterile transitions). This last case is 
especially well-motivated theoretically, since it constitutes one of the simplest ways 
to reconcile i the ANA with other puzzles in the neutrino sector such as the solar 
neutrino problem as well as the LSND result i and the possible need for a few eV 
mass neutrino as the hot dark matter in the Universe i. 

The main aim of this talk is to compare the and the —>■ z/^ transitions 

using the the new sample corresponding to 535 days of the Super-Kamiokande data. 

*To Appear in Proceeding of the RINGBERG EUROGONFERENGE - NEW TRENDS IN NEU¬ 
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This analysis uses the latest improved calculations of the atmospheric neutrino fluxes 
as a function of zenith angle, including the muon polarization effect and taking into 
account a variable neutrino production point & 


2. Atmospheric Neutrino Oscillation Probabilities 

The expected neutrino event number both in the absence and the presence of 
oscillations can be written as: 
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and Paf 3 is the oscillation probability of —>■ for given values of E,^, cosOy and h, 
i.e., Pa 0 = Pi^a Ey, cos9y, h). In the case of no oscillations, the only non-zero 

elements are the diagonal ones, i.e. Paa = 1 for all a. 

Here rit is the number of targets, T is the experiment’s running time, Ey is the 
neutrino energy and is the flux of atmospheric neutrinos of type a = /r, e; Ep is 
the final charged lepton energy and e{Ep) is the detection efficiency for such charged 
lepton; a is the neutrino-nucleon interaction cross section, and 6y is the angle between 
the vertical direction and the incoming neutrinos (cos 0,^=1 corresponds to the down¬ 
coming neutrinos). In Eq. (^, h is the slant distance from the production point to 
the sea level for a-type neutrinos with energy Ey and zenith angle 6y. Finally, Ka is 
the slant distance distribution which is normalized to one i. 

The neutrino fluxes, in particular in the sub-GeV range, depend on the solar 
activity. In order to take this fact into account in Eq. (I). a linear combination of 
atmospheric neutrino fluxes and which correspond to the most active Sun 
(solar maximum) and quiet Sun (solar minimum) respectively, is used. 

For definiteness we assume a two-flavor oscillation scenario, in which the os¬ 
cillates into another flavour either z/^ —^ z/g , r'/i —>■ z^s or z/^ ^ z/^. The Schrodinger 
evolution equation of the z/^ — ux (where X = e, r or s sterile) system in the matter 
background for neutrinos is given by 
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Here Gp is the Fermi constant, p is the matter density at the Earth, M is the nncleon 
mass, and Ye {Yn) is the electron (nentron) fraction. We dehne Am^ = — m\ in 

snch a way that if Am^ > 0(Am^ < 0) the nentrino with largest mnon-like component 
is heavier (lighter) than the one with largest X-like component. For anti-nentrinos 
the signs of potentials Vx shonld be reversed. We have nsed the approximate analytic 
expression for the matter density prohle in the Earth obtained in ref. 0. In order to 
obtain the oscillation probabilities Pap we have made a nnmerical integration of the 
evolntion eqnation. The probabilities for nentrinos and anti-nentrinos are different 
becanse the reversal of sign of matter potential. Notice that for the case there 

is no matter effect while for the Ug case we have two possibilities depending 

on the sign of Am^. For Am^ > 0 the matter efects enhance neutrino oscillations 
while depress antineutrino oscillations, whereas for the other sign (Am^ < 0) the 
opposite holds. The same occnrs also for z/^ —z/g. Althongh in the latter case one 
can also have two possible signs, we have chosen the most nsnally assnmed case where 
the mnon nentrino is heavier than the electron nentrino, as it is theoretically more 
appealing. Notice also that, as seen later, the allowed region for this sign is larger 
than for the opposite, giving the most conservative scenario when comparing with 
the present limits from CHOOZ. 


3. Atmospheric Neutrino Data Fits 

Here I describe onr £t method to determine the atmospheric oscillation parameters 
for the varions possible oscillation channels, including matter effects for both —>■ Ve 
and —>■ Vg channels. The steps required in order to generate the allowed regions of 
oscillation parameters were given in ref. 0. I will comment only that when combining 
the results of the experiments we do not make use of the double ratio, 
but instead we treat the e and /i-like data separately, taking into account carefully 
the correlation of errors. It is well-known that the double ratio is not well suited from 
a statistical point of view due to its non-Gaussian character. Thus, following ref. 00 
we define the as 

x" = YkK'" - <“"■*) ■ - Np"’), (4) 
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where I and J stand for any combination of the experimental data set and event-type 
considered, i.e, / = (A, a) and J = {B, (3) where. A, B stands for Frejus, Kamiokande 
sub-GeV, IMB,... and a,P = e, p. In Eq. (|D is the predicted number of events 
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Fig. 1. Xmin fixed Am^ versus Am^ for each oscillation channel for Super-Kamiokande sub-GeV 
and multi-GeV data, and for the combined sample. Since the minimum is always obtained close to 
maximum mixing the curves for i/g for both signs of coincide. 

calculated from Eq. (2) whereas is the number of observed events. In Eq. (^) 

(Jtheory error matrices containing the experimental and theoretical 

errors respectively. They can be written as 

= ^a{A) papiA, B) a/3 (5), (5) 

where pQ,p{A,B) stands for the correlation between the a-like events in the A-type 
experiment and /3-hke events in S-type experiment, whereas (Ja{A) and (^/^{B) are the 
errors for the number of a and /3-hke events in A and B experiments, respectively. 

We compute Pa^{A, B) as in ref. 0. A detailed discussion of the errors and corre¬ 
lations used in the analysis can be found in Ref.0. We have conservatively ascribed a 
30% uncertainty to the absolute neutrino flux, in order to generously account for the 
spread of predictions in different neutrino flux calculations. Next we minimize the 
function in Eq. and determine the allowed region in the sin^ 29 — /Aw? plane, for 
a given conhdence level, dehned as, 

X' = xLn + 4.61 (9.21) for 90 (99)% C.L. (6) 

In Fig. 1^ we plot the minimum (minimized with respect to sin^ 26) as a function 
of Aw?. Notice that for large Am^ ^0.1 eV^, the is nearly constant. This happens 
because in this limit the contribution of the matter potential in Eq (^) can be neglected 
with respect to the Aw? term, so that the matter effect disappears and moreover, 
the oscillation effect is averaged out. In fact one can see that in this range we obtain 
nearly the same y^ for the —>■ Vr and —*■ Vg cases. For very small Am^^l0“'^ eV^, 
the situation is opposite, namely the matter term dominates and we obtain a better £t 
for the z/^ —>■ channel, as can be seen by comparing the —>■ Vj. curve of the Super- 
Kamiokande sub-GeV data (dotted curve in the left panel of Fig. |l]) with the 
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and —*• Vf. curves in the left panel of Fig. [^). For extremely small Am^ ^ 10“^ 
eV^, values is quite large and approaches a constant, independent of oscillation 
channel, as in the no-oscillation case. Since the average energy of Super-Kamiokande 
multi-GeV data is higher than the sub-GeV one, we hnd that the limiting Am^ value 
below which approaches a constant is higher, as seen in the middle panel. Finally, 
the right panel in Fig. |I| is obtained by combining sub and multi-GeV data. A last 

Table 1. Minimum value of and the best fit point for each oscillation channel and for different 
data sets. For > Vg the minimum is practically independent of the sign of Am^ as the 
minimum is located at maximum mixing angle. 


Experiment 




t'/. ^ 

Super-Kam 

A.min 

7.1 

8.2 

7.3 

sub-GeV 

Am^ ( 10“^eV^ ) 

0.11 

1.9 

1.2 


sin^ 26 

1.0 

1.0 

0.97 

Super-Kam 

Y^ ■ 

A.min 

6.3 

7.9 

10.8 

multi-GeV 

Am? ( 10“^eV^ ) 

1.5 

3.5 

24.7 


sin^ 26 

0.97 

1.0 

0.72 

Super-Kam 

Y^ ■ 

'^min 

14.3 

16.8 

21.8 

Gombined 

Am^ ( 10“^eV^ ) 

1.6 

2.6 

1.5 


sin^ 26 

1.0 

1.0 

0.97 

All experiments 

Y^ ■ 

A-min 

47.2 

48.6 

48.6 

Gombined 

Am^ ( 10“^eV^ ) 

2.9 

3.5 

3.0 


sin^ 26 

1.0 

1.0 

0.99 


point worth commenting is that for the —»• case in the sub-GeV sample there 
are two almost degenerate values of Am^ for which attains a minimum while for 
the multi-GeV case there is just one minimum at 1.5 x 10“^eV^. Finally in the third 
panel in Fig. ^ we can see that by combining the Super-Kamiokande sub-GeV and 
multi-GeV data we have a unique minimum at 1.6 x 10“^eV^. 

4. Results for the Oscillation Parameters 

The results of our £t of the Super-Kamiokande sub-GeV and multi-GeV at¬ 
mospheric neutrino data are given in Fig. In this figure we give the allowed region 
of oscillation parameters at 90 and 99 % GL. One can notice that the matter effects 
are similar for the upper right and lower right panels because matter effects enhance 
the oscillations for neutrinos in both cases. In contrast, in the case of > z/^ with 
Am^ < 0 the enhancement occurs only for anti-neutrinos while in this case the effect 
of matter suppresses the conversion in z/^’s. Since the yield of atmospheric neutrinos 
is bigger than that of anti-neutrinos, clearly the matter effect suppresses the overall 
conversion probability. Therefore we need in this case a larger value of the vacuum 
mixing angle, as can be seen by comparing the left and right lower panels in Fig. p. 














Fig. 2. Allowed regions of oscillation parameters for Super-Kamiokande for the different oscillation 
channels as labeled in the figure. In each panel, we show the allowed regions for the sub-GeV data 
at 90 (thick solid line) and 99 % CL (thin solid line) and the multi-GeV data at 90 (dashed line) 
and 99 % CL (dot-dashed line). 


Notice that in all channels where matter effects play a role the range of acceptable 
Am^ is shifted towards larger values, when compared with the —>■ Vr case. This 
follows from looking at the relation between mixing in vacuo and in matter. In fact, 
away from the resonance region, independently of the sign of the matter potential, 
there is a suppression of the mixing inside the Earth. As a result, there is a lower cut 
in the allowed Am^ value, and it lies higher than what is obtained in the data fit for 
the —>■ i^r channel. 

It is also interesting to analyse the effect of combining the Super-Kamiokande 
sub-GeV and multi-GeV atmospheric neutrino data. Gomparing the results obtained 
with 535 days given in the table above with those obtained with 325 days of Super- 
Kamiokandea we see that the allowed region is relatively stable with respect to the 
increased statistics. However, in contrast to the case for 325.8 days, now the z/^ —>• i^r 
channel is as good as the —>■ when only the sub-GeV sample is included, with a 
clear Super-Kamiokande preference for the z/^ —> Uj. channel. As before, the combined 
sub-GeV and multi-GeV data prefers the z/^ —>• ux, where X = r or sterile, over the 
solution. 

To conclude this section I now turn to the predicted zenith angle distributions 
for the various oscillation channels. As an example we take the case of the Super- 
Kamiokande experiment and compare separately the sub-GeV and multi-GeV data 
with what is predicted in the case of no-oscillation (thick solid histogram) and in all 















Figure 3: Angular distribution for Super-Kamiokande electron-like and muon- like 
sub-GeV and multi-GeV events together with our prediction in the absence of oscil¬ 
lation (dot-dashed) as well as the prediction for the best fit point for —>• Vg (solid 
line), —>• z/e (dashed line) and > z/,- (dotted line) channels. The error displayed 
in the experimental points is only statistical. 


oscillation channels for the corresponding best fit points obtained for the combined 
sub and multi-GeV data analysis performed above (all other histograms). This is 
shown in Fig. ^ 

It is worthwhile to see why the z/^ —>■ z/g channel is bad for the Super-Kamiokande 
multi-GeV data by looking at the upper right panel in Fig. |^. Glearly the zenith 
distribution predicted in the no oscillation case is symmetrical in the zenith angle 
very much in disagreement with the data. In the presence of z/^ —>■ z/g oscillations 
the asymmetry in the distribution is much smaller than in the z/^ —>■ 1 /^ or z/^ —>■ Ug 
channels, as seen from the figure. Also since the best fit point for z/^ —> Ug occurs at 
sin(26*) = 1, the corresponding distributions are independent of the sign of Am^. 

5. Atmospheric versus Accelerator and Reactor Experiments 

I now turn to the comparison of the information obtained from the analysis of the 
atmospheric neutrino data presented above with the results from reactor and acceler¬ 
ator experiments as well as the sensitivities of future experiments. For this purpose I 
present the results obtained by combining all the experimental atmospheric neutrino 
data from various experiments^. In Fig. §we show the combined information obtained 








































from our analysis of all atmospheric neutrino data involving vertex-contained events 
and compare it with the constraints from reactor experiments such as Krasnoyarsk, 
Bugey, and CHOOzi, and the accelerator experiments such as CDHSW, CHORUS, 
and NOMAD 0. We also include in the same hgure the sensitivities that should be 
attained at the future long-baseline experiments now under discussion. 

The hrst important point is that from the upper-right panel of Fig. §one sees that 
the CHOOZ reactoiil data already exclude completely the allowed region for the 
Vf. channel when all experiments are combined at 90% CL. The situation is different 
if only the combined sub-GeV and multi-GeV Super-Kamiokande are included. In 
such a case the region obtained is not completely excluded by GHOOZ at 90% GL. 
Present accelerator experiments are not very sensitive to low Am^ due to their short 
baseline. As a result, for all channels other than —>■ v^. the present limits on neutrino 

oscillation parameters from GDHSW, GHORUS and NOMAD 0 are fully consistent 
with the region indicated by the atmospheric neutrino analysis. Future long baseline 
(LBL) experiments have been advocated as a way to independently check the ANA. 
Using different tests such long-baseline experiments now planned at KEK (K2K) 0, 
Fermilab (MINOS) 0 and GERN ( IGARUS 0, NOE 0 and OPERA 0) would test 


the pattern of neutrino oscillations well beyond the reach of present experiments. 
These tests are the following: r appearance searches, NCjCC ratio which measures 

, and the muon disappearance or CCnear/CCfar test. The second test 

can potentially discriminate between the active and sterile channels, i.e. —>■ Vj. 

and —>■ However it cannot discriminate between —>■ z/^ and the no-oscillation 
hypothesis. In contrast, the last test can probe the oscillation hypothesis itself. Notice 
that the sensitivity curves corresponding to the disappearance test labelled as KEK- 
SK Disappearance at the lower panels of Fig. ^ are the same for the z/^ —> Uj. and the 
sterile channel since the average energy of KEK-SK is too low to produce a tau-lepton 
in the far detector. In contrast the MINOS experiment has a higher average initial 
neutrino energy and it can see the tau’s. Although in this case the exclusion curves 
corresponding to the disappearance test are in principle different for the different 
oscillation channels, in practice, however, the sensitivity plot is dominated by the 
systematic error. As a result discriminating between z/^ —>■ z/^ and z/^ —z/^ would be 
unlikely with the Disappearance test. 

In summary we hnd that the regions of oscillation parameters obtained from the 
analysis of the atmospheric neutrino data on vertex-contained events cannot be fully 
tested by the LBL experiments, when the Super-Kamiokande data are included in 
the ht for the z/^ —> channel as can be seen clearly from the upper-left panel of 
Fig. § One might expect that, due to the upward shift of the Am^ indicated by 
the ht for the sterile case, it would be possible to completely cover the corresponding 
region of oscillation parameters. This is the case for the MINOS disappearance test. 
But in general since only the disappearance test can discriminate against the no¬ 
oscillation hypothesis, and this test is intrinsically weaker due to systematics, we hnd 
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Fig. 4. Allowed oscillation parameters for all experiments combined at 90 (thick solid line) and 
99 % CL (thin solid line) for each oscillation channel as labeled in the figure. We also display the 
expected sensitivity of the present accelerator and reactor experiments as well as to future long- 
baseline experiments in each channel. The best fit point is marked with a star. 






















that also for the sterile case most of the LBL experiments can not completely probe 
the region of oscillation parameters indicated by the atmospheric neutrino analysis. 
This is so irrespective of the sign of Am?\ the lower-left panel in Fig. shows the 
^ channel with Am? < 0 while the —>■ case with Am? > 0 is shown in the 

lower-right panel. 
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